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We systematically develop the full quantum theory for the electromagnetic induced transparency 
(EIT) and slow light properties in ultracold Bose and Fermi gases. It shows a very different property 
from the classical theory which assumes frozen atomic motion. For example, the speed of light in- 
side the atomic gases can be changed dramatically near the Bose-Einstein condensation temperature, 
while the presence of the Fermi sea can destroy the EIT effect even at zero temperature. From exper- 
imental point of view, such quantum EIT property is mostly manifested in the counter-propagating 
excitation schemes in either the low-lying Rydberg transition with a narrow line width or in the 
D2 transitions with a very weak coupling field. We further investigate the interaction effects on 
the EIT for a weakly interacting Bose-Einstein condensate, showing an inhomogeneous broadening 
of the EIT profile and nontrivial change of the light speed due to the quantum many-body effects 
beyond mean field energy shifts. 

PACS numbers: 42.50.Gy, 67.85.-d, 32.80.Ee 



I. INTRODUCTION 

Electromagnetic induced transparency (EIT) has been 
extensively studied in quantum optics [lj, |2| and applied 
to serve as quantum memory elements through light- 
atom interactions 0] . It is also realized in ultracold Bose 
gas 3 to slow down the light speed in an unprecedented 
way. The transparency along with significant group ve- 
locity delay opens up great opportunities of light storage 
and retrieval in neutral atomic systems, as one of the 
candidates in quantum information processing. 

In the classical theory of EIT, it usually assumes frozen 
atomic motion which is valid only in limited conditions 
when the recoil momenta from laser interactions are neg- 
ligible. When the thermal cloud at several hundreds fiK 
or higher is considered, the frequency shift resulted from 
the atomic motion is comparable to the EIT transparency 
window. The Doppler effect can cause decoherence to the 
probe field due to the large spread of atomic velocity 
distribution Q. In the other regime of ultracold atoms 
considered in this article, a Gross-Pitaevskii (GP) equa- 
tion has been used in Bose-Einstein condensates (BEC) 
in the zero temperature limit to investigate light propa- 
gation dynamics @ . The dispersion relations of light are 
unaffected due to the extremely low energy atomic dy- 
namics (< kHz) even in the counter-propagating scheme, 
providing the great opportunity and experimental flexi- 
bility in the quantum information application based on 
the EIT effect. However, this is not the case for exam- 
ple when Rydberg states are considered in the applica- 
tion Q , in which the recoil energy of light imparting on 
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atoms is comparable to the inverse of lifetime for Rydberg 
transitions. Also at finite temperature across the quan- 
tum degeneracy, atomic Bose or Fermi statistics starts 
to play a crucial role in determining EIT property where 
GP equation is not applicable. 

To develop the general quantum theory of EIT, in this 
paper we include the atomic kinetic energy explicitly to 
get the dark state dynamics. Compared to the energy 
scale of the coupling field and the decay rate of the ex- 
cited state, the recoil kinetic energy is found to be sig- 
nificant mostly in the counter-propagating scheme either 
in the low-lying Rydberg transition or the D2 transition 
of a very weak field. For bosonic gases, the group veloc- 
ity of the transmitted light is strongly modified near the 
BEC transition temperature. For fermionic gases, on the 
other hand, there is no EIT allowed due to the presence 
of a Fermi sea. We also calculate the EIT profile for co- 
propagating scheme, and find almost no effects of such 
quantum statistics, and EIT is robust in a wide range of 
temperature (even across the quantum degenerate tem- 
perature) for both bosonic and fermionic atoms. Finally, 
we also include interaction effects in bosonic gases. We 
find that, besides of the simple mean-field energy shift, 
the EIT profile can be broadened by the many-body ef- 
fects, significantly increase the group velocity when the 
interaction is stronger. Our results show that the EIT ex- 
periment can be a good candidate to explore the quantum 
many-body properties of the underlying quantum gases. 

In this article, we review the classical theory on EIT 
in Sec. II and derive the general quantum theory that 
includes atomic dynamics and its quantum statistics in 
alkali atomic gases in Sec. III. We demonstrate the EIT 
property in noninteracting Bose and Fermi gases using 
ultraviolet transitions of low-lying Rydberg states and 
D2 transitions for both co- and counter-propagating ex- 
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FIG. 1: (Color online) Electromagnetic induced transparency 
(EIT) in A-type atoms, (a) The atoms are interacting with 
two laser fields either through co (dashed)- or counter (solid)- 
propagating scheme, (b) Control Q.2 and probe fields fii cou- 
ple two hyperfine ground states \a) and \b) with the excited 
state |e) with detunings A2 and Ai respectively. F is the spon- 
taneous emission rate of |e). (c) Linear susceptibility x times 
i versus probe light detuning Ai in unit of T. Absorption 
(solid-blue) and dispersion (dashed-red) profiles are plotted 
from the classical theory of EIT in resonance of the control 
field A 2 = 0. 



citation schemes. In Sec. IV, we apply the mean field 
approach to the EIT in a weakly interacting BEC and 
conclude in Sec. V. 



II. CLASSICAL THEORY ON EIT 

The classical theory assumes frozen atomic motion so 
that the physics of EIT mainly relies on the quantum in- 
terference of two laser fields coupled with A type atomic 
structure as shown in Fig. [TJ The atoms interact with 
control and probe fields for transitions \b) — > |e) and 
I a) — > |e) of Rabi frequencies Q2 and &i respectively, ei- 
ther in co- or counter-propagating excitation scheme. Sin- 
gle photon detunings are defined as Ai = lo\ — (u e — co a ), 
A 2 = L02 — (<-*J e ~ w b)- The relatively weak probe field 
fii propagates into a medium interacting with a strong 
control field ^2, which modifies the energy of the ex- 
cited state via AC Stark shift. The one dimensional wave 
equation for slow- varying probe electric field E(z,t) can 
be described by [|[ 



up a self-consistent Maxwell-Bloch equation which fully 
accounts for atomic coherence and light propagation dy- 
namics. Use P(z,t) = (pd a e&ae(z,t)) with the atomic 
density p, dipole matrix element d ae , and slow- varying 
atomic coherence operator a ae (z,t) in a rotating frame 
of probe light central frequency uj\ and wavevector k\, 
we may express susceptibility and group velocity in res- 
onance of the constant control field (A 2 = 0) as [1| 
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where g = d ae u>\ / (2heoV) is field coupling constant 
with the quantization volume V. iV| f7| 2 associates with 
optical depth (OD) that determines the reduction of 
group velocity and N is the total number of atoms. The 
group velocity is derived from v g = duji/dk(uji) where 
fc(wi) = wi/c + Im [ifcix(wi)/(2eo)]. In the above, we 
specifically show the results at the zero absorption pint 
where Ai = 0. The absorption and dispersion profiles 
are shown in Fig. 1(c) for the Rabi frequency VL2 — T. 
The typical results from classical theory show that the 
width of transparency window is determined by 25^2 with 
a normal dispersion in it. The group velocity at zero ab- 
sorption point depends on the control field strength and 
optical depth, which decreases as 1^2 decreases or with a 
larger OD. Another useful property of oscillator strength, 
will be used later to characterize 
the coupling strength between S and P level transitions 
of frequency u> and dipole matrix element d, and m e is 
the mass of electron with charge e. 

The above results delineate the essence of induced 
transparency through electromagnetic couplings from 
classical theory. In the next section, we generalize the 
treatment of EIT to include the atomic dynamics and 
quantum statistics especially in nonintcracting quantum 
degenerate gases at finite temperature. 
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where erj is the vacuum permittivity and the atomic po- 
larization in Fourier space is P(z,uj) — x( z > UJ )E(z, uj) 
with the electric susceptibility x(z,u;). Along with the 
optical Bloch equations of atoms 8], the above makes 



III. NONINTERACTING QUANTUM 
DEGENERATE GASES AT FINITE 
TEMPERATURE 



In general, the Hamiltonian (H = H a+H al) of nonin- 
teracting quantum degenerate gases in momentum space 



3 



is 



^ = E^-( a P«p+^p+^pj + 

p p 
p 

Hal = -hn 1 ^2el +pi a p -hn 2 ^2el +p2 b p + h.c. 7 (3) 



where Ha includes the kinetic and internal state energy 
for atoms, and Hal represents the atom-light interac- 
tions. Dipole approximation, — d • E, and rotating wave 
approximation (RWA) have been made to the atom-light 
interactions, pi 2 are recoil momenta of the probe and 
control fields respectively. Note that we have used inter- 
action picture for atom and light fields to remove the free 
evolution terms (we keep the atomic kinetic energy), and 
atom fields satisfy the bosonic or fermionic commutation 



relations, e.g. [a p , dL]± = 6 ] 



We may express the Hamiltonian in terms of the basis, 
A p = [d p , b p+Prl e p+Pl ] , with p, = p x p 2 , and let 
H = J2 P Ap ' M p ' A p , where the matrix M p is (h = 1) 
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In general, we can diagonalize if by a similarity trans- 
formation S p such that in the new basis B p = S p A p , 
it becomes a diagonal matrix A p = S P M P S~ . We are 
interested in the regime when the probe field is weak so 
we can solve the above Hamiltonian in the leading or- 
der of small Hi. Under the adiabatic approximation, the 
atoms initially prepared on the ground state |o) follow 
these light fields and form a dark state, $ P which is im- 
mune to spontaneous emission of the excited state (up to 
a normalization constant), 
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where ezj(p) and e B± {p) are the eigenvalues of the un- 
perturbed Hamiltonian for the dark and two bright states 
respectively, 
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where A 2 = A 2 + (p + p r ) 2 / (2m) and 
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ground and excited states coupled by the control field. 

This dark state provides the ingredient to the disper- 
sion and absorption of light propagation in the atoms. 
For a nonintcracting quantum degenerate gas, we derive 
the linear susceptibility x( w i) under the dark state basis 
assuming a constant control field, 
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The angle <fi p represents the mixing of the hyperfinc where 



$(r) is the atomic field operator in real space and is 
the direction perpendicular to pi with the cross section 
area A. The expectation value of Eq. flSJ) is calculated 
from the dark state in Eq. ([3]). Note that x( w i) has spa- 
tial dependence of z which is implicit because we assume 
a uniform atomic density p. The boson/fermion number 
distribution is np p — l/[exp{[e£>(p) — /x]/(fcsT)} =p 1] 
where ks is the Boltzmann constant, and the chemi- 
cal potential p, is determined by the number conserva- 
tion X] p n /3 P = N. For convenience we define a ker- 
nel function, F p which will be encountered quite often 
throughout this paper. The second line of F p is derived 
by adding a phenomenonological spontaneous decay rate 
(r) of the excited state into the noninteracting Hamil- 
tonian. The dispersion relation from the above quantum 
theoretical results indicates a frequency shift for single 
and two-photon detunings from atomic central (p) and 
recoil (pi and p r ) momentum distributions. The absorp- 
tion width of probe field is mainly determined by the 
energy scale of natural decay rate T. Note that Eq. © 
reduces to the classical results when the atomic motion is 
assumed frozen. In subsequent subsections, we study non- 
interacting bosonic/fcrmionic gases and investigate the 
EIT property when the excited states are low-lying Ry- 
dberg state and D2 transitions for both co- and counter- 
propagating excitation schemes. 



Results of counterpropagating excitation 
scheme 



From Eq. ©, there are three energy scales, 0,2, I\ 
and E r = pf/(2m), that play the major roles in the dis- 
persion and absorption profiles of EIT. When Dl or D2 
transition of the alkali atoms is considered, the sponta- 
neous decay rate (r) is in the range of MHz which exceeds 
low energy atomic dynamics (< kHz). The quantum the- 
ory of EIT in the previous section reduces to the clas- 
sical results as if the atomic motion is frozen when the 
control field strength SI2 is also in the order of T. How- 
ever, the spontaneous emission rate of the Rydberg state 
can be in the range of kHz H and even below for larger 
principal quantum number [7]. This motivates us to use 
Rydberg transitions to investigate the quantum effects 
of EIT where O2 ~ E r rj T. For comparison, we also 
calculate D2 transitions when using a weak control field, 
SI2 ~ E r <C r, which can show quite different properties 
compared to classical results. 



1. Low-lying Rydberg bosons 

Here we first consider the counterpropagating excita- 
tion scheme as shown in Fig. QJa) an d choose the excited 
levels of rubidium (Rb) and potassium (K) as low-lying 
Rydberg states. The radiative lifetimes for low-lying Ry- 
dberg atoms (nP levels) [9j can reach to the comparable 
recoil energy scale when the principal quantum number 
is around n — 24 [lfj. Note that the dipole-dipole inter- 
action of Rydberg atoms can be significant and relevant 
here. To simplify the discussion in this work, we only 
consider the situation that the dipole-dipole interaction 
or dipole blockade is weak, e.g. low-lying Rydberg ex- 
citations (n < 30) [llj]- The long-range van der Waals 
interaction coefficient | Ce | for 30P-30P rubidium asymp- 
totes is in the range of (0.26 - 4.7) x 10 16 a.u. pi that 
gives us the blockade radius R c = 0.88 — 1.29 [l3l.ll4T| 
for a density p = 10 14 cm~ 3 with the effective excitation 
Rabi frequency SI = 2ir x 100 kHz. This radius provides 
an upper estimation of the strength of dipole blockade. 
For potassium atoms (n = 30), we have blockade radius 
R c = 1.3 — 1.67 /im if the same atomic density and exci- 
tations are applied. The dipole-dipole interaction energy 
can be further reduced by decreasing SI so that it is not 
relevant when low-lying Rydberg states are concerned. In 
the future, we can also include the dipole-dipole interac- 
tion of Rydberg atoms. 

For a noninteracting Bose gas, the transition temper- 
ature can be found analytically for a three dimensional 
uniform gas T c = 27rft 2 (p/C(3/2))( 2 / 3 ) /(mk B ) \'i where 
£ is Riemann zeta function, and p is atomic density of 
mass to. The chemical potential approaches the mini- 
mum of eigenenergy ec(p = 0) = Ai when temperature 
approaches T c . We find T c w 0.4 pK for 87 Rb of a density 
p = 10 14 cm- 3 . 

At T < T C) the condensate particle density (p c ) is 
derived from subtraction of the excited particle density 
p c (T) — p — Pex- The electric susceptibility of finite tem- 
perature are resulted from the condensate and non con- 
densate parts, 
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where the first term is exactly the classical result in Eq. 
@ if Pc — Pi and the second term is thermal depletion. 
To calculate the group velocity for the probe light, we 
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FIG. 2: (Color online) Counterpropagating excitation scheme 
for EIT in bosons using the Rydberg transition (I 1 " 1 = 26^is) 
of 87 Rb. The control field is SI2 = ir with the detuning 
A2 = which is used throughout the article if not specified. 
Absorption (Re[i%]) (solid-blue) and dispersion (Im[ix]) (dash- 
red) profiles are plotted at different temperatures (a) T=5, 
(b) 0.5, (c) 0.35, and (d) 0.1/iK where transition temperature 
T c ~ 0.4pK (straight dotted line in (e)). Dash-dotted black 
line guides the eye to the zero of the plots. In (e), the group ve- 
locity v g is plotted for various control fields, ^2/^ = 0.5(O) > 
!(□), 2(0), 4(+), 10(») over temperatures in log scale. The 
detuning is set as A2 = — p^/(2m) so v g is calculated at the 
transparency condition of Ai = 0. 



need the oscillator strengths (/) between the ground and 
low-lying Rydberg states, |21P) of K and |24P) of Rb, 
which are 3.1 x 1CT 7 and 2.8 x 1CT 6 respectively [jj]. The 
branch ratio of transitions to P1/2 are large en oug h that 
P3/2/-P1/2 ~ 4 and 8 for K and Rb respectively [rTlliljj so 
we may take the values of oscillator strengths as a correct 
order of magnitude. For D2 transitions of potassium and 
rubidium atoms, / = 0.658 and 0.682, much larger than 
the low-lying Rydberg states. 

In Fig. [5J we show the electric susceptibility (x) and 
group velocity (v g ) of the probe field in 87 Rb. The excited 
state is chosen as low- lying Rydberg states |e) = |24P 3 / 2 ) 
with spontaneous decay time = 28.3 fj,s. When 

T > T c in (a) and (b), there is no transparency due 
to the large bosonic momentum distribution that devi- 
ates from the resonance condition of free absorption. The 
relaxation due to the Doppler width at high tempera- 
ture can be estimated as kBTh\2/h. Just below T c in 
(c), two absorption (Re[ix]) lines appear, and the trans- 
parency condition starts to emerge at Ai w 3r that is 
due to the recoil energy shift (E r ) for the probe field 
wavelength 299 nm. When T is much smaller than T c 
in (d), the condensate component dominates over the 
thermal one, and we have a perfect EIT with zero ab- 
sorption. The two absorption peaks can be identified as 
Ai = {£1 + E r ± [{Ex - E r f + 4h 2 nl}^ 2 }/(2h) where 
Ei = pf/(2m). Note that the absorption peak in plot 
(a) when T >• T c can be identified as Ai = E\j% 
due to the thermally-averaged transparency condition 
(Ax = E r /h). 

In Fig. HJe), the group velocity is plotted for vari- 
ous control fields from 0.5 to 10r. For small we find 



FIG. 3: (Color online) Counterpropagating excitation scheme 
for EIT in bosons using D2 transition (r _1 = 26ns) of 87 Rb. 
Absorption (Re[ix])( s °hd-blue) and dispersion (Im[ix]) (dash- 
red) profiles are plotted as a comparison to Fig. [2] at the tem- 
peratures (a) T=5, (b) 0.5, (c) 0.35, and (d) 0.1/xK. Dash- 
dotted black line guides the eye to the zero of the plots, 
and Q.2 = O.Oir, A2 = 0. In (e), the group velocity v g 
is plotted for fi 2 / r = 0.003(solid-blue), O.Ol(dashed-green), 
0.03(dash-dotted red) over temperatures in log scale across 
T c (straight dotted line). 



a significant reduction of group velocity when T < T c , 
and it saturates in the limit of T — > 0. There is no slow 
light at a higher temperature because no transparency is 
supported if O2 is too small. This contrasts to the con- 
ventional experiments using Dl or D2 transitions of Rb 
where f2 2 is large 0,13. From Fig. He), we demonstrate 
the interplay between the control fields VL2 and temper- 
ature T for the effect of EIT on slow light propagation. 
We find that this may be an efficient and precise method 
to determine the T c in bosonic quantum gases. 



2. D2 transition bosons 

In Fig. [3l we demonstrate the similar results compared 
to Fig. [5]when the D2 transition (r _1 = 26 ns) is chosen 
with a very weak control field ^2 = O.Oir in (a)-(d). The 
sharp transparency window on top of the large absorption 
width appears when temperature is lower, which signifies 
a reduction of momentum distribution. In (e), the group 
velocity behaves similarly to Fig. ^e) at the low T limit 
which saturates as atomic motions do not matter any- 
more, while at higher T we still see slow light up to the 
temperature of mK (not shown) though EIT is very weak 
due to the weak control field we choose here. The rela- 
tively large T of D2 transition makes the decoherence of 
Doppler effects possible only when T > HT/kB- 



3. Fermions 

For a noninteracting fermionic gas, it is straightforward 
to derive the EIT property by using fermionic number 
distribution in Eq. ©. In Fig. [H we show the elec- 
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FIG. 4: (Color online) Counterpropagating excitation scheme 
for EIT in fermions. (a,b) Absorption (solid-blue) and disper- 
sion (dash-red) profiles are plotted for 40 K at temperatures 
of 5 and 0.5 (iK (Fermi temperature Tf = 1.98/xK). (c) Two 
energy bands of free and shifted-recoil particles are shown 
as p 2 /(2m) and — Ai + (p + p r ) 2 /(2m) respectively (Ai is 
chosen to be equal to recoil energy for example). The en- 
ergy difference between two energy bands is demonstrated for 
fermions where solid-red arrows represent different excitation 
paths from the Fermi surface, (d) Same energy bands of (c) 
are shown for bosons. The zero momentum (p = 0) state (red 
filled circle) match these two bands exactly for Bose-Einstein 
condensation. 

trie susceptibility of probe field in 40 K atoms. The ex- 
cited state is chosen as low- lying Rydberg states |21P 3 / 2 ) 
(transition wavelength 288 nm) with spontaneous decay 
time 25.4 [is. In Fig. |4](a) and (b), we find no EIT across 
the Fermi temperature Tf = 1.98 /iK which indicates a 
mismatch of resonance condition for the probe field as 
shown in (c). This mismatch comes from the Fermi dis- 
tribution of atoms that has a Fermi surface of momentum 
and the recoil energy shift, E r . It contrasts with bosonic 
atoms that mostly distributes around zero momentum as 
shown in (d) when condensation sets in. The absorption 
peaks in (a) and (b) can also be identified by the re- 
coil energy of probe laser, Ai = Ex/%. The effect of the 
temperature is the absorption width that is larger for a 
higher temperature. Similar to Rydberg bosons, we may 
increase f2 2 to overcome the momentum mismatch that 
EIT emerges. For D2 transitions, we have similar results 
as bosons in Fig. [3](a,b) where the transparency window 
shows up on top of the large absorption profile when a 
weaker control field is used (f2 2 <§C T). As the control field 
increases, e.g. f2 2 ~ R the EIT profile does not change 
across Tf as if atomic motions are frozen. 

B. Results of copropagating excitation scheme 

Here we consider the copropagating excitation scheme 
where the recoil momentum is reduced to the energy 
transfer of hyperfine splitting of the ground states which 
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FIG. 5: (Color online) Copropagating excitation scheme for 
EIT in bosons. Absorption (solid-blue) and dispersion (dash- 
red) profiles of 87 Rb are shown in the same temperature con- 
ditions (a)-(d) as in Fig. [5] In (e), the group velocity v g is 
plotted for various control fields, Q, 2 /V = 0.5(O), !(□), 2(0) 
over temperatures from 0.01 to 1 0/iK in log scale, which shows 
significant reduction and saturates in the low T limit. High 
temperature limit has more reduction due to the narrower 
transmission window. 

are 6.834 and 1.285 GHz for 87 Rb and 40 K respectively. 
The main difference between co- and counter-propagating 
excitation schemes is the magnitude of recoil energy {E r ) 
that decides the transparency condition. From Eqs. © 
and ([9]), we can see that when p r 0, F p is almost inde- 
pendent of p, so that the temperature dependence of the 
EIT effects will not be expected significantly. 

First we consider low-lying Rydberg bosons and study 
the temperature dependence of EIT in Fig. [5j When 
T < T c as shown in Fig. [5jc) and (d), it is similar to the 
counterpropagating scheme where the probe field inter- 
acts largely with condensation particles that have zero 
momentum distribution. There are again two absorption 
peaks similarly identified as Rydberg bosons in the coun- 
terpropagating excitation scheme with E r ss 0, and the 
transparency condition is Ai « due to negligible E r . 
When T ^> T c , the absorption widths are broadened due 
to atomic motions, which makes the transmission window 
narrower as in (a) . In (e) we calculate the group velocity 
at the transparency condition. It is plotted over a range 
of temperature from 0.01 to 10 in log scale for vari- 
ous control fields from 0.5 to 2r. It saturates in the low 
T limit, the same as the counterpropagating scheme, but 
becomes smaller at higher T due to the narrower trans- 
mission window. 

In Fig. [6l we consider low-lying Rydberg fermions, and 
investigate the EIT property and group velocity. Absorp- 
tion peaks can be identified in the same way as bosons 
in copropagating scheme. For lower temperature in (b), 
the widths of absorption peaks are narrower indicating 
a momentum distribution near the Fermi surface. They 
can be estimated in the order of 2^/EpEi/H from Eq. 
(|9]) where Ep is the Fermi energy. The separation be- 
tween the peaks approaches 2fi 2 when the control field 
is stronger. The transmission window becomes narrower 
when T is higher, which decreases the group velocity as 
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FIG. 6: (Color online) Copropagating excitation scheme 
for EIT in fermions. Absorption (solid-blue) and dispersion 
(dash-red) profiles of 40 K. Same temperature conditions are 
used as in Fig. 2] Group velocity v a is plotted for tempera- 
tures from 0.01 to IO^iK in log scale in (c) for various control 
fields, n 2 /F = 0.5(O), !(□)> 2(0). It saturates in the low T 
limit and more reduction of velocity is observed due to the 
narrower transmission window at higher T. 



shown in Fig. [6jc). Contrary to the counterpropagating 
scheme, there is always a transparency window for the 
temperatures considered here with f2 2 = T. It is due to 
negligible recoil energy shift that two-photon resonance 
condition is reached even for large atomic momentum 
distribution at high temperature. 

For D2 transitions with a larger T, similar to the 
Rydberg bosons/fermions, there is a transparency win- 
dow due to negligible recoil energy shift for f2 2 down to 
O.OOir. And again, the EIT appears on top of the absorp- 
tion profile as in the counterpropagating scheme with D2 
transitions. 



IV. WEAKLY INTERACTING BOSE-EINSTEIN 
CONDENSATE 



where in the first expression we show in general all the 
possible atom-atom interactions U st — 4irh 2 a st /rn with 
s-wave scattering lengths a s t in our A type atomic config- 
uration (s, t = a,b or c). In the second expression we use 
the Hartree-Fock approximation |15l | for the interactions 
of U a b, U ae , and p c = (ao) 2 /V is the mean field conden- 
sate density. The neglected terms include Ubb, U ee , Ube, 
and other interactions between the atomic fields b p , e p , 
and a Pi AQ. They are at least 0(£l\) which are negligible 
in the weak field limit. We then absorb the mean field 
energy shifts of the bare atomic fields b p and e p into the 
nonintcracting Hamiltonian of Eq. ([3]), which becomes 
(cf. Eq. @D) 



M' p = 


= PcUael 4 





A j_ (P + Pr) 



2m + PcUab 

z 2m 



-o 2 



A 2 





(P+Pr) 



2m 



-Oo 



-n a 
(p+pi) 2 

2m 



PcU a 

, (12) 



where we have extracted the common mean field energy 
p c U ae that we derive the effective noninteracting Hamil- 
tonian with shifted detunings Ai = Ai — pJJ a e and 
A 2 = A 2 + p c (Uab — U ae ) in the last line. 

Using the dark state (3 p introduced in Sec. 3 and with 



Eq. ([12]). we let a p « j3 p and 



fNnSr: 



/Wo 



where (3 p ^o are quasiparticles. We then apply Hartree- 
Fock-Bogoliubov and Popov approximations for the in- 
teractions of U aa in Eq. (TTT1) . and diagonalize the to- 
tal Hamiltonian by Bogoliubov transformation. Define 
Hd = Ha + Hal + Hjj and introduce the chemical po- 
tential p to take care of the total number conservation, 
we have 



Here we include the interaction energy in the noninter- 
acting Hamiltonian of Bose gas and investigate the EIT 
in a weakly interacting Bose-Einstein condensate (BEC). 
The interaction Hamiltonian is 



Hu = ^ [^^P+q^P'-q^P'^P 



p,p',q 

Ubb~bi:~bl,b„ r b 



bfc°P + qV-qV°P t" ^eeep+q'v-qVep 
i J2 [ 2C/ ^P+q« P '-q«P'^P 

pp'q 

2t r aee P +q 6 P'-qVe P + 2C4 e e p+q & P ,_ q Ve p 

a p+q«l'-qV«P + PcU ab ^ b l b P 



pp'.q 



+ p c Uae^e. p e p , 



(11) 



H D - pN = E g - /iAT c + ^e(q)a q a q , (13) 
where 



U aa , 



Eg = e D (p = 0)N C + -^(i\H + 2Ni x ) + 



, N c NU ae 



2V 



V 



q>0 



(14) 



is the ground state energy including the mean field energy 
and energy correction, and cd(p) = Ai + p 2 /(2m). The 
eigenvalue and eigenstates of the Hamiltonian are e(q) = 
(ef (q) — e 2 )^ 2 and a q = cosh6>(q)/3 q +sinh6>(q)/3L q where 
tanh26»(q) = e 2 /ei(q), ei(q) = eo(q) - p + 2p c U aa , and 
(2 = PcUaa- The chemical potential is p — hAi + p c U aa 
which is found by regarding gapless Bogoliubov excita- 
tion energy, e(q = 0) = 0. 

Using the obtained new ground state wavefunction, we 
derive the electric susceptibility for the probe field in the 
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weakly interacting BEC, 



Xc(wi) + XQd(ui) + Xth{ui), 
P 



+ — ^>6Cosh26>] ; 



(15) 



p#0 



where we separate the condensate(C), quantum de- 
pletion (QD), and thermal depletion (TH) contribu- 
tions to The shifted kernel function is F p with 
Ax — > Ax and A 2 — > A 2 . Bose distribution is npjq) = 
l/[ e e (q)/C £ BT) _ j^j j? or Bogoliubov quasi-particles [15j . 

Without loss of generality, we let U a b — U aa , and as- 
sume U ae ~ (which is just to have an excitation en- 
ergy shift). In Fig. [7J we demonstrate the EIT prop- 
erty for counter- and co-propagating excitation fields 
in weakly interacting Bose-Einstein condensates using 
low-lying Rydberg transitions. We choose the scatter- 
ing length a s seven times of 106ao for 87 Rb where do 
is Bohr radius, and explore the effects of quantum de- 
pletion (p ex R3 O.lp with negligible thermal depletion at 
T = lOnK) on the profiles. The overall profiles in (a) and 
(c) bear the main features for noninteracting Bose gas 
below T c in Fig. HJd) and EJd) respectively. The mean 
field interaction energy, p c U aa ~ 0.9r, shifts the trans- 
parency conditions of noninteracting BEC, which we de- 
note as the double arrows in (a) and (c). It is of kHz 
that can be observable in experiments. Two absorption 
peaks are still resolved and the transparency conditions 
are mainly determined by the condensate particles. The 
insets from the Bogoliubov excitations can be seen simi- 
larly as in Fig. [2{b) andHfb) respectively at temperature 
above T c . XQd(w\) and xth(w\) are contributions from 
finite momenta which broaden the shapes of absorption 
peaks. This interaction-induced broadening in interact- 
ing Bose-Einstein condensates can be observed through 
increasing the interaction energy via the Feshbach res- 
onance. In Fig. [JJb) and (d), we calculate the group 
velocity near the transparency point and find that it de- 
creases (increases) for co(counter)-propagating excitation 
scheme when the interaction strength increases. It is due 
to the narrower (broader) transmission window from the 
Bogoliubov particle contributions. For even stronger in- 
teraction strengths, the Bogoliubov approximation fails. 



V. CONCLUSION 

We theoretically investigate the electromagnetic in- 
duced transparency in ultracold atomic gases. The 
atomic dynamics and quantum statistics significantly 
modify the light propagating properties in quantum 
degenerate gases concerning noninteracting bosons, 
fermions, and weakly interacting Bose-Einstein conden- 
sate for co- and counter-propagating excitation schemes. 
Throughout the systematic studies of EIT property, we 
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FIG. 7: (Color online) EIT property in weakly interacting 
S7 Rb condensate. Dispersion (dash-red) and absorption (solid- 
blue) profiles are shown in (a) counter- and (c) copropagating 
excitation schemes. The respective group velocities (b) and 
(d) are calculated at the transparency conditions and plot- 
ted over interaction strengths (a s ). The scattering length and 
temperature are chosen as a s — 7 x 106 ao and T — lOnK in 
(a) and (c) where the insets are i(xqd + Xth) from Bogoli- 
ubov particles p ex ~ O.lp. The coupling atomic levels and the 
atomic density p are the same as in Fig. O and the double 
arrows in (a) and (c) are the mean field interaction energy 

Shifts, pcUaa- 



carefully analyze the interplay of three energy scales: 
control field strength, temperature, and spontaneous 
emission rate. We use low-lying Rydberg excited 
states and D2 transitions to demonstrate the EIT prop- 
erty across the characteristic temperatures (T c ,Tf) of 
noninteracting bosons and fermions. In the counter- 
propagating scheme, we show a transition of the EIT 
property across T c of bosons when the control field is 
close to the recoil energy. It can be used as a sensitive 
determination of T c . For fermions at low temperature, 
the presence of the Fermi sea can destroy the EIT ef- 
fect, and the absorption widths can be determined by the 
Fermi energy. As the control field increases, the EIT pro- 
file can be described by the classical theory as if atomic 
motions are frozen. The results of D2 transitions are sim- 
ilar to those of the Rydberg states. In the copropagating 
scheme, EIT is robust for both bosons and fermions in a 
wide range of temperature even with a weak control field 
due to negligible recoil energy. 

We further apply the Hartree-Fock-Bogoliubov mean 
field approach to the weakly interacting Bose-Einstein 
condensate, and find that EIT is modified by quantum 
and thermal depletions due to excitation particles. The 
mean field energy from atom-atom interactions shifts 
the transparency condition, and an interaction induced 
broadening is shown in absorption widths which can be 
observable by tuning atomic scattering lengths via Fesh- 
bach resonance. The group velocity is also calculated and 
analyzed in details for various schemes of different tem- 
peratures, control field or interaction strengths. Future 
developments might involve using EIT to investigate the 
many-body and strongly interacting atomic system, and 
it may provide a fruitful information of atomic dynam- 
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ics or synthesize magnetic gauge fields to simulate 
spin-orbit coupling quantum degenerate gases [20| . 
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